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ABSTRACT 

We have analysed all the good quality XMM-Newton data publicly available for the bright 
ULXs HolmberglX X-l and NGC 1313 X-l, with the aim of searching for discrete emission 
or absorption features in the Fe K band that could provide observational evidence for the 
massive outflows predicted if these sources are accreting at substantially super-Eddington 
rates. We do not find statistically compelling evidence for any atomic lines, and the limits 
that are obtained have interesting consequences. Any features in the immediate Fe K energy 
band (6-7 keV) must have equivalent widths weaker than ^30 eV for HolmberglX X-l, and 
weaker than ^50 eV for NGC 1313 X-l (at 99 per cent confidence). In comparison to the sub- 
Eddington outflows observed in GRS 1915+105, which imprint iron absorption features with 
equivalent widths of ~30 eV, the limits obtained here appear quite stringent, particularly when 
Holmberg IX X- 1 and NGC 1 3 1 3 X- 1 must be expelling at least 5-10 times as much material 
if they host black holes of similar masses. The difficulty in reconciling these observational 
limits with the presence of strong line-of-sight outflows suggests that either these sources are 
not launching such outflows, or that they must be directed away from our viewing angle. 
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1 INTRODUCTION 

Ultraluminous X-ray Sources (ULXs) are extra-nuclear point 
sources in external galaxies observed to be more luminous in 
X-rays than the Eddington luminosity Le for a stellar mass 
(~10Mq) black hole (Lx > 10' i9 erg s _1 ); on rare occasions 
ULX luminosities have been observed to exceed 10 41 erg s _1 



dWalton et all 1201 lbl : ISutton et ail l2012h an d even to reach as 
high as ~ 10 erg s" - 1 i lFarrelletalJl2009t) . This combination 
has led to extended debate over the nature of these sources. The 
majority of early theories were based around either the pres- 
ence of intermediate mass blac k hole s (IMBHs: 10 2 < A/ BH < 
10 5 Mm : IColbert & Mushotzkvll 19991) . supe r-Eddington accretion 
states for standard mass binary systems (e.g. IPoutanen et al]|2007l . 
iFinke & B6ttchej|2007h or the inferre d luminosities be ing artifi- 
cially high due to anisotropic emission jKing et all200lh . although 
observations of excited emission lines, typically He II A4686, 
from either the accretion disc of the ULX or its surrounding neb- 
ula appear to rule out highly beamed emission via photon co unt- 
ing arguments (e.g.lPakull & MirionHl2002l . rKaaret et alj|2004 and 



iBerghea et all20H 

More recently, proposed interpretations for bright ULXs 
(Lx ~ 10 40 erg s" 1 ) tend to incorporate elements from each of 
these three ideas, invoking black holes only slightly larger than 
those seen in Galactic binary systems (Mbh < 100 Mq) with 

* E-mail: dwalton@ast.cam.ac.uk 



mildly super-Eddington accretion rates (luminosities up to a few 
times Le), such that the accretion disc is geometrically thicker 
than t he standard thin discs pred icted for more moderate accretion 
rates (shakura & Sunvaevl 1973) and forces so me mild anisotropy 
('slim' discs; see e.g. lAbramowicz et alj|l980l) . For recent reviews 
on the observa tion al status and the pote ntial nature of ULXs see 
iRobertsj J2007h and lFeng & Sorij feoill) . 

One of the observational results in recent X-ray studies of 
ULXs that has sparked significant interest is that many of these 
sources display curvature in their X-ray continuum above ~3 keV 
dStobbart et afl|2006t iGladstone et~ai1l2009l ; IWalton et al.ll201 lab . 
Although some of these are lower luminosity sources that dis- 
play (hot) thermal disc-like spectra, in which high energy curva- 
ture would naturally be expected, even the brighter ULXs with two 
apparently distinct emission components, potentially analogous to 
the disc and Comptonised corona observed in Galactic black hole 
binaries (BHBs), also appear to show curvature in their high en- 
ergy components. The hard, Comptonised emission rarely displays 
similar curvature in the stand ard sub-Eddington accretion sta tes of 
BHBs (for a recent review see lRemillard & McClintockl 20061) . and 
if this is an intrinsic difference in the high energy continuum, it 
could imply that there are fundamental differences between the ac- 
cretion onto ULXs and their lower luminosity BHB cousins. Given 
that the basic accretion geometry is not expected to evolve with 
black hole mass, this might indirectly suggest we are not viewing 
sub-Eddington accretion. 

In this high-Eddington framework, IGladstone et al. I <2009h 
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Table 1. Basic details of the two bright ULX s co nsidered in t his work, HolmberglX X-l and NGC 1313 X-l, and the XMM-Newton observations analysed. 
Distances are taken from lPaturel et alJfeOOA and lTullvl lll988) respectively. 



Source 


RA 


DEC 


D 


z 


OBSID 


Obs. Date 


Duration 


Good Exposure 




(h:m:s) 


(d:m:s) 


(Mpc) 








(ks) 


(EPIC-pn; ks) 


HolmberglX X-l 


09:57:53.2 


+69:03:48.3 


3.55 


0.000153 


0112521001 


10/04/2002 


1 1 


7 












0112521101 


16/04/2002 


12 


8 












0200980101 


26/09/2004 


119 


88 


NGC 1313 X-l 


03:18:20.0 


-66:29:11.0 


3.7 


0.001568 


0106860101 


17/10/2000 


42 


27 












0150280601 


08/01/2004 


55 


10 












0150181101 


16/01/2004 


9 


6 












0205230201 


01/05/2004 


13 


3 












0205230301 


05/06/2004 


12 


9 












0205230601 


07/02/2005 


14 


10 












0405090101 


15/10/2006 


123 


99 



propose that the high energy curvature could be due to Comptonisa- 
tion in an optically thick corona, which shrouds the inner accretion 
disc. The effect of the corona causes the observed temperature of 
the disc (identified here with the soft component) to appear artifi- 
cially low as only the outer disc is observed directly. Low disc tem- 
peratures have frequently been inferred for bright ULXs, which if 
otherwise associat ed with the inner d isc would imply the presence 
of IMBHs (see e.g. iMiller et al.l2004 ). This interpretation has been 
slightly modified bv lMiddleton et al" J201 ll) . who propose that the 
cool, quasi-thermal emission previously associated with the disc is 
instead thermalised emission from a photosphere associated with 
the base of an outflowing wind, as strong outflows are ubiquitously 
predicted b y models for high and partic ularly super-Eddington ac- 
cretion (see lOhsuga & Mineshigdfeoill . and references therein). 

Galactic BHBs frequently display evidence for mass outflows, 
which at moderately high accretion rates (during the thermal dom- 
inated states) take the form of equatorial disc winds with out- 
flow velocities v out < 1000 km s -1 . Prominent examples are the 
BHBs GRO J1655- 40 jMiller et al.ll2006t iDfaz Trigo et alj|2007t) 
and G RS 1915+105 iKotani et al.l2000l : lLee et al.l2002klUeda etaU 
Hill Neilsen & Le d l2009t) . When these outflows impinge upon 
our line of sight to the central source, their observational conse- 
quence is to imprint absorption features onto the intrinsic X-ray 
continuum, typically in the form of narrow absorption lines from 
highly ionised species, with the most prominent features arising 
from the Ka transitions of highly ionised iron (Fe XXV and/or 
XXVI). In the case of GRS 191 5+105, perhaps the best studied BHB 
outflow. lNeilsen & Led J2009I) show that the strength of the absorp- 
tion features generally increases as the luminosity of the source, 
and hence its Eddington ratio increases. The mass of this source i s 
relatively well known (M B h = 14 ± 4 Mm : iGreiner et al.ll200lh . 
and it is observed to ra diate up to, and possibl y slightly in excess 
of its Eddington limit JVierdavanti et alj|2010h . While the geom- 
etry of the BHB disc winds observed at moderate accretion rates 
are largely equatorial, the scale height of these wind s is generally 
expected to incre ase with Eddington ratio (see e.g. lAbramowiczl 
l2005llKindl200l) . so at higher accretion rates the outflows should 
probably cover a larger solid angle. 

If instead, super-Eddington outflows are present in ULXs but, 
despite their potential large scale heights, do not occur along our 
line of sight, rather than viewing absorption features we should in- 
stead observe reprocessed emission from the material in the out- 
flow. Indeed, reprocessed emission lines consistent with neutral 
(or moderately ionised) iron are seen ubiquitously in the spectra 



of Ga lactic high mass X-ray binaries (HMXBs; see lTorreion et al .1 
l20ld) . which are expected to accrete largely from the stellar winds 
launched by their massive companions. Furthermore, the majority 
of ULXs themselves are widely expected to be HMXB analogues, 
owing to the their a pparent correlation with recent star formation 
dSwartz et alj|2009h . Recent population studies of ULXs appear 
to support this picture, as their luminosity function appears con- 
sistent with being a smooth extrapolation of the HMXB luminos- 
ity function (at least for sources in spiral/irregular type g alaxies ; 
iGrimm et al j2003llWalton et alj|201 lbllSwartz et ai1l201 ll) . 

Here, we present an analysis of the Fe K band for two bright 
(Lx ~ 10 40 erg s _1 ) ULXs with hard X-ray spectra, Holmberg IX 
X-l (also known as M81 X-9) and NGC 1313 X-l, in which we 
search for evidence of atomic iron features that could be associated 
with outflowing material. These sources have been selected as they 
represent the datasets with the highest quality data at high energies 
for such bright, isolated ULXs. The paper is structured as follows: 
section[2]describes our data reduction prodecure, and section[3]de- 
scribes the analysis performed. We discuss our results in section [4] 
and summarise our conclusions in section[5] 



2 OBSERVATIONS AND DATA REDUCTION 

In this work we make use of all the pu blically available XM M- 
Newton obs ervations with good EPIC-pn dStriider et al .Il200lh and 
EPIC-MOS l lTurner et al J 200 ll spectra of the two bright ULXs con- 
sidered, HolmberglX X-l and NGC 1313 X-l. In the latter case, 
there are a number of observations pointed at NGC 1313 X-2 in 
which X- 1 fell on a chip gap in one or more of the XMM-Newton 
detectors which were not included in our analysis. In addition, we 
only consider observations in which the respective sources are de- 
tected reliably above background up to 10 keV. Data reduction 
was carried out with the Science Analysis System (SAS vl 1.0.0) 
largely according to the standard prescription provided in the on- 
line guidefl 

Observation data files were processed using EPCHAIN and 
EMCHAIN to produce calibrated event lists for the EPIC-pn and 
EPIC-MOS detectors respectively. Spectra were produced sepa- 
rately for each observation by selecting only single and double 
events (single to quadruple events) for EPIC-pn (EPIC-MOS) us- 
ing XMMSELECT, and periods of high background were treated ac- 

1 http://xmm.esac.esa.int/ 
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Table 2. Key parameters obtained for HolmberglX X-l and NGC 1313 X-l with the phenomenological ultraluminous state continuum model (see text). 
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Source 


Obs. 


N H 
(10 21 cm" 2 ) 


kT in 
(keV) 


kT c 
(keV) 


T 


X 2 (x 2 /d-o.f.) 


HolmberglX X-l 


Combined 
0200980101 


1.5 ±0.1 
1.5 ±0.1 


0.23 ± 0.01 
0.23 ±0.01 


2.4 ±0.1 
2.3 ±0.1 


8.7 ±0.3 
9.2 ±0.4 


1.08 (1998/1851) 
1.02(1804/1762) 


NGC1313X-1 


Combined 
0405090101 


2.5 ±0.1 

2 fi+ ' 2 
z -°-o.i 


0.22 ± 0.01 
0.22 ± 0.01 


9 q+0.2 
" a '°-0.1 
9 9+0-1 
z - z -0.2 


7.6 ±0.4 
8.5 ±0.5 


1.16(1867/1608) 
1.08 (1415/1305) 



cording to the method outlined bv lPiconcelli et alj d2004) with the 
signal-to-noise ratio maximised for the full energy band considered 
(0.3-10.0 keV). Circular source regions were used, and although 
the exact size varied from observation to observation depending on 
the location of chip gaps, etc. these were typically of radius ~30- 
40". Larger circular regions typically of radius ~50" and ~75" free 
of other sources and on the same CCD as the ULX were selected 
in order to sample the background for the EPIC-pn and EPIC- 
MOS detectors respectively. Redistribution matrices and auxiliary 
response files were generated with RMFGEN and ARFGEN. After 
performing the data reduction separately for each of the observa- 
tions, the spectra were then combined using the FTOOI0 ADDSPEC 
in order to maximise the data quality in the Fe K band. Finally, 
EPIC-pn and EPIC-MOS spectra were re-binned using GRPPHA to 
have a minimum of 25 counts in each energy bin, so that the prob- 
ability distribution of counts within each bin can be considered 
Gaussian, and hence the use of the \ 2 statistic is appropriate when 
performing spectral fits. 

The observations that met our selection criteria and are hence 
included in our analysis are listed in Table [TJ the resulting good 
exposure times obtained in the stacked spectra are 103 and 238 ks 
with the EPIC-pn and combined EPIC-MOS detectors respectively 
for HolmberglX X-l, and 162 and 381 ks for NGC 1313 X-l, al- 
though these are largely dominated by the single full orbit observa- 
tions available for each source (contributing 88/200 ks to the EPIC- 
pn/ EPIC-MOS spectra of HolmberglX X-l, and 90/215 ks to the 
spectra of NGC 1313 X-l respectively). 



model the soft comp onent and cool, o ptically thick thermal Comp- 
tonisation (COMPTTi lTitarchukl 19941) to model the curvature of the 
hard component. The inner radius of the disc is fixed at 6 Rg (= 
GM/c 2 ) and the seed photon temperature for COMPTT is linked to 
the inner temperature of the disc. Both components are modified 
by ne utral absorption with solar abundances (TBABSi lwilms et"al] 
l200d) . The overall quality of fit obtained over the full 0.3-10.0 keV 
band is generally fairly good, with xl = 1998/1851 and 1867/1608 
for the stacked spectra of HolmberglX X-l and NGC 1313 X-l re- 
spectively, but if we focus on the quality of fit obtained with these 
models (without refitting) at higher energies (2-10 keV) we find 
it to be much improved, with xl = 1429/1394 and 1218/1151. We 
therefore conclude that the continuum is well modelled over the en- 
ergy band of interest for this study (~5-9 keV) for both sources, 
see Fig.HJ The results obtained for the key continuum parameters 
are quoted in Table|2] For comparison, we also repeat our analysis 
for the long observations of each source separately; we obtain good 
agreem ent with the results initially presented in iGladstone et al.1 
d2009h . 

To confirm whether it is reasonable to combine the individual 
observations into a single spectrum, we also apply this continuum 
model to each of the observations individually. For each source we 
model all the individual observations considered simultaneously, 
requiring that the neutral absorption remains the same for each ob- 
servation, but allowing the other key parameters to vary (see Ap- 
pendix|A]for full details). The key parameters obtained for each of 
the observations are all quite similar, therefore we conclude it is 
reasonable to combine the observations included in this work. 



3 SPECTRAL ANALYSIS 

Throughout this work, spectr al modelling is performed with 
XSPEC vl2.7.0r dArnau j U9961) . and all quoted uncertainties are 
the 90 per cent confidence limits for a single parameter of inter- 
est, unless stated otherwise. In the following, HolmberglX X-l 
andNGC 1313 X-l are modelled over the full 0.3-10.0 keV energy 
range. The EPIC-pn and EPIC-MOS spectra are modelled simulta- 
neously with all parameters tied between the two spectra, although 
we do include a multiplicative constant in the EPIC-MOS spectrum 
in order to account for possible cross-calibration uncertainties be- 
tween the detectors. This value is always found to be within ~10 
per cent of unity. 

3.1 Continuum Model 

We begin by modelling the stacked spectra of these sources with a 
simple phenomenological interpretation of the ultraluminous st ate 
concept, including a cool disc (DISKPN ; iGierlihski et aljH~999l) to 

2 http://heasarc.nasa.gov/ftools/ftools_menu.html 



3.2 Fe K Band 

Having successfully reproduced the shape of the continuum for 
these two sources, we now proceed to search for evidence of narrow 
iron emission or absorption features which might indicate a large, 
super-Eddington outflow. We are primarily interested in iron in this 
work owing to its high cosmic abundance and fluorescent yield, and 
also owing to the fact that it is difficult to fully ionise. Furthermore, 
the XMM-Newton response function varies slowly and smoothly 
over the iron Ka bandpass. To search for these features, we add 
a narrow (a — 10 eV) Gaussian feature with a normalisation that 
can either be positive to model potential emission features, or nega- 
tive to model potential absorption features. The restframe energy of 
this component (taking into account the redshift adopted for each 
source, as quoted in TableQJ is systematically varied across the en- 
ergy range 5-9 keV in 100 steps, in order to incorporate the key 
iron transitions and also (generously) allow for any potential blue 
shifts, as would be expected for absorption features from a strong 
outflow. This a pproach is broadly sim ilar to that adopted in the re- 
cent works by iTombesi et al.l J2010al fbh during their initial search 
for potential emission and absorption features in the high energy 
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Energy (keV) Energy (keV) 



Figure 1. Continuum model fits to the stacked spectra of Holmberg IX X-l (left panel) and NGC 1313 X-l (right panel). The models consist of a cool disc 
co mponent (red) and cool , optically thick Comptonisation (blue), in accordance with the ultraluminous state concept for super-Eddington accretion proposed 
bv lGladstone et all j2009h . The lower panels show the data/model ratio for the models constructed over the energy range of interest for this work, 5-9 keV; 
these models clearly provide a good representation of the high energy continuum. 



spectra of active galaxies. At each energy increment we record the 
Ax 2 improvement provided by the inclusion of the Gaussian fea- 
ture, as well as the equivalent width (EW) of any potential feature 
at that energy, and its 90 and 99 per cent confidence limit These 
are obtained with the EQWIDTH command in XSPEC by generating 
10000 parameter sets based on the best fit model and the covariance 
matrix, which contains information on the parameter uncertainties 
obtained with this fit, and extracting the 90 and 99 per cent uncer- 
tainty limits from the distribution of equivalent width values ob- 
tained from each of these parameter sets. 

The results obtained for Holmberg IX X-l and NGC 1313 X- 
1 are shown in Fig. [2] and Fig. [3] respectively; in each case the 
left panel displays the results obtained with the stacked spectrum, 
and the right panel displays the results obtained when considering 
only the longest observation of the source. The top panels display 
the A\ 2 improvement as a function of line energy, multiplied by 
the sign of the normalisation to differentiate between emission and 
absorption, and the lower panels display the limits on the equiv- 
alent widths obtained. For clarity, we highlight the energies of the 
Ka transitions of neutral iron, iron XXV and iron XXVI with dotted 
vertical lines. We also show EW — ±30 eV with dashed horizon- 
tal lines, which roughly repres ent the strongest absorp tion features 
observed in GRS 1915+105 bv lNeilsen & Lee] j2009l) , although as 
noted by those authors the EW values obtained in that work are 
actually likely to be lower limits o wing to their contin uum fitting 
technique, and in GRO J1655-40 bv lMiller et ai]d2006h . 

First of all, we wish to stress that we have not found any sta- 

3 Note that these are confidence limits on line strength, which are not 
strictly the same as the upper limits on the strengt h of lines that shoul d 
be detected at a given detection confidence level (seejKashyap et al1 l2010l) , 
but are simpler to calculate in this instance, and should most likely provide 
more conservative estimates for the maximum strength of any lines intrin- 
sically present. 



tistically significant detections of any features, either in absorption 
or in emission. Although a few energies are picked out by our anal- 
ysis as offering the most improvement, this is never greater than 
A\ 2 ~ 7 over the energy range considered. Therefore our dis- 
cussion is restricted to the EW limits that can be placed instead, 
although these are interesting to consider in their own right. We 
find that the data are of sufficient quality to exclude very strong 
(EW > 50 eV) iron features in either absorption or emission 
across a large portion of energy range immediately around the Fe K 
transitions, although the constraints naturally become less restric- 
tive at very high energies (>7. 5-8.0 keV) as the data quality de- 
clines. The high energy data quality is better for Holmberg IX X-l, 
so the limits that can be placed are strongest in this case, where any 
persistent atomic features in the spectrum below ~7.5 keV must 
have an equivalent width less than ~30 eV. In the case of NGC 1313 
X-l, a moderately strong (EW < 50 eV) emission feature consis- 
tent with neutral iron and absorption feature consistent with iron 
XXVI are both still permissable. In the following section we dis- 
cuss the implications of these limits on the potential presence of 
large outflows and the accretion geometry for these ULXs. 



4 DISCUSSION 

We have analysed the available XMM-Newton observations for 
two bright ULXs with hard X-ray spectra, Holmberg IX X-l and 
NGC 1313 X-l, with the purpose of searching for direct observa- 
tional evidence of massive super-Eddington outflows via atomic 
features in the iron K band. In order to maximise the high energy 
data quality we stacked the available data with simultaneous EPIC- 
pn and EPIC-MOS coverage, and proceeded to search for narrow 
features in either absorption or emission. We do not find any statis- 
tically compelling evidence for any narrow features, so instead we 
examine the limits that can be placed on any features intrinsically 
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5678956789 
Restframe Energy (keV) Restframe Energy (keV) 

Figure 2. Top panels: the A\ 2 improvement provided by including a narrow (a = 10 eV) Gaussian line, as a function of line energy, for Holmberg IX X-l. 
The normalisation of the line is allowed to be either positive or negative to allow for both emission and absorption features, and the Ax 2 improvement is 
multiplied by the sign of the best fitting line normalisation to indicate whether the improvement is provided by an emission or an absorption feature. We find 
no statistically compelling evidence for any narrow iron features at high energies. Bottom panels: 90 (blue) and 99 (red) per cent confidence contours for the 
equivalent width of the included narrow line as a function of line energy, indicating the line strengths any persistent narrow features intrinsically present could 
have and still remain undetected with the current data. The rest frame transitions of neutral, helium-like and hydrogen-like iron (6.4, 6.67 and 6.97 keV) are 
shown with vertical dotted lines, and for comparison we also plot dashed horizontal lines representing EW = ± 30 eV. This analysis is repeated for both the 
stacked spectrum (left panels) and the longest single XMM-Newton observation (right panels). 




5678956789 
Restframe Energy (keV) Restframe Energy (keV) 

Figure 3. As in Fig.|2] but here for NGC 1313 X-l. 
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present in the spectrum that remain undetected. For comparison 
we also analyse the full orbit observations available for these two 
sources separately, which make up ~85 and 55 per cent of the total 
good exposure time analysed for Holmberg IX X-l and NGC 1313 
X-l respectively. In the latter case the additional observations con- 
tribute more significantly to the stacked spectrum, and a more not- 
icable improvement is obtained with their inclusion. However, the 
results are similar in both cases, and we focus our discussion on the 
limits obtained with the stacked spectra. 

4.1 Absorption Lines 

First we consider the limits that can be placed on persistent ab- 
sorption features, given the strong outflows expected from super- 
Eddington accretion. For Holmberg IX X-l any absorption features 
in the immediate Fe K band (6-7 keV) must have an equivalent 
width less than ~30 eV, based on the more conservative 99 per cent 
limits obtained. The constraints are not so strong for NGC 1313 
X-l, where an absorption feature consistent with Fe XXVI could 
potentially still have an equivalent width of up to ~50eV. The 
implication is that any line-of-sight outflows present cannot have 
any greater impact on the spectrum of Holmberg IX X-l than 
the strongest sub-Eddington Galactic BHB outflows, observed in 
GRS 1915+105 and GRO J1655-40, and can only have at most a 
moderately stronger impact for NGC 1313 X-l. Furthermore, we 
also draw attention to the outflow in the Galactic BHB candidate 
IGRJ17091-3624 dicing et alj|2012l). This source dis plays strong 
similarities to GRS 1915+105 jAltamirano et al.l201lT) . and though 
its distance and mass are not currently known, through this compar- 
ison it has tentatively been suggested that it could also be emitting 
at close to, or even in excess of its Eddington limit. Although the 
statistical significances of the iron features in this source are not 
particularly high, if real their reported equivalent widths are ~40 
and ~90eV for Fe XXV and XXVI when using a cool, optically 
thick Comptonisation model for the continuum (similar to that ob- 
served in GRS 1915+105, and to the ultraluminous state interpreta- 
tion). Any undetected features in the immedite Fe K band must be 
much weaker than these. 

However, assuming black hole masses of 10 Mq for Holm- 
berg IX X-l and NGC 1313 X-l, these sources are observed to per- 
sistently radiate at ~5—10 times their Eddington limits. In this sce- 
nario, we would expect them to be ejecting at least 5-10 times 
the amount of material as a simila r source accreting at the Ed- 
dington limit (i.e. GRS 1915+105). IPonti et alj J2012h find that, 
even under the assumption of a constant wind opening angle, the 
ratio of the mass outflow rate to the innermost accretion rate, 
M on t/M a CC , increases with increasing Eddington ratio, so increas- 
ing the luminosity by a factor of 10 increases the amount of ma- 
terial expelled by an even larger factor (assuming this increase 
in luminosity relates to an equivalent increase in A/ acc ). Given 
that EW scales roughly linearly with absorbing column density, 
if the ionisation state and abundances of the winds remained con- 
stant we might therefore have expected to see iron absorption fea- 
tures with EW > 300 eV if these sources are at ~10 times the 
Eddington limit and their associated outflows are along our line- 
of-sight, assuming solar abundances. For t he Z ~ 0.5 Zq abun- 
dances of the host galaxies in these cases llMakarova etalll2002l ; 
lHadfield & Crowtherll2007l : IPintore & Zampierill2012h , the expec- 
tation would still be EW > 150 eV. 

If these outflows are truly absent, this would argue very 
strongly against a highly super-Eddington interpretation for these 
ULXs. Both sources display hard X-ray spectra with evidence for 



two separate continuum emission components, similar to some of 
the complex 'intermediate' sub-Eddington states, with very cool 
disc temperatures (~0.2 keV) obtained when associating the soft 
component with the accretion disc. However, both sources also dis- 
play evidence for the high energy curvature in their hard compo- 
nents that might distinguish ULXs from the typical sub-Eddington 
states seen in BHBs, and is potentialy t he key signature of th e 
proposed 'ultraluminous' accretion state dGladstone et alj|2009h . 
In fact, GRS 1915+105 may show some similar spectral charac- 
teristics, including the high energy curvature, w hen radiating close 
to its Eddington limit l lVierdavanti et al j|2010h . although the pres- 
ence of a reflection component in this source complicates the is- 
sue. If the predicted outflows for strong super-Eddington accretion 
really are absent, this curvature might still indicate another high 
Eddington-ratio accretion state, but one that is still roughly Ed- 
dington limited rather than strongly super-Eddington, and observed 
when L/Le ~ 1. This would suggest the black hole masses of 
Holmberg IX X-l and NGC 1313 X-l are ~50-100 M . 

There are, however, some important caveats to these limits that 
must be considered, especially given the ionised nebulae that many 
bright ULXs are observed to reside within. In some cases, includ- 
ing Holmberg IX X- 1 , the outer regions of these nebulae appear to 
be more consistent with being ionised by shocks instead of X-ray 
photoionisation, and these regions are observed to be expandin g 
JPakull & MMorj2002fcrAbolmasov etai]| 20071 : Icleh et alj|2012l) . 
Where present, these shock-ionised regions have generally been in- 
terpreted as potential evidence for the large scale impact on the 
ISM of strong winds launched from the central ULX (although, of 
course, the mere presence of an outflow from the accretion disc 
does not by itself imply super-Eddington accretion). Historically, 
interpretations in which they are instead the expanding remnants of 
the explosive even t in which the ULX was formed have also been 
proposed (see e.g. iRoberts et al.ll2003h . but giv en the typical age s 
of the surrounding stellar populations (see e.g. iGrise et al]|201ll ). 
the nebulae appear more consistent with the wind inflated scenario. 
Therefore, we now consider whether, despite the absorption limits 
obtained from our X-ray analysis, we could still be viewing these 
bright ULXs through strong, super-Eddington outflows. 

First, we note briefly that when considering the stacked spec- 
tra we are strictly obtaining limits on persistent features, and that 
the outflows observed in BHBs do vary with time, broadly corre- 
lating with accretion state ( see e.g. Miller et alJl200fj . lMiller et al .1 
120081 iNeilsen & Leeil200l. IPonti et alj|2012h . However, although 
the sources do display some mild long term flux variability be- 
tween the observations considered here, it is certainly not clear that 
these variations are related to accretion state changes, particularly 
given that the individual spectra obtained from the observations in- 
cluded in this work all appear intrinsically similar (see appendix lAl). 
Therefore, a search for persistent features seems reasonable. Fur- 
thermore, similar (if mildly weaker) constraints are obtained with 
the single full-orbit observations of e ach source, in w hich virtually 
no short term variability is observed jHeil et alj2009h . 

Although the limits obtained across the immediate Fe K band 
(6-7 keV) are fairly stringent, they naturally become less restric- 
tive at higher energies. Stronger absorption features could still be 
present in the spectrum were they to arise in winds with very 
high outflow velocities. This is particularly true for NGC 1313 
X-l where persistent narrow absorption features with equivalent 
widths up to ~75eV (or even ~100eV at truly extreme veloci- 
ties) could still currently be present. In contrast, even high velocity 
features are constrained to have EW < 50 eV for Holmberg IX 
X-l. In either case, outflow velocities of « ou t > 20,000 kms~ 



© 2010 RAS, MNRAS OOO.mfiTI 



Searching for Winds in ULXs: the Fe K band 7 



(0.08c) would be required in order to invoke absorption features 
significantly stronger than currently possible in the immediate Fe 
K band. To date, the fastest potential outflow observed in any 
BHB has a velocity o f v out ~ 10,000 km s" 1 (IGRJ17091- 
3624; lKingetal.ll2012h . although much more typical velocities 
are v ou t < 1000 kms -1 . However, there is growing evidence 
for w i nds with n out > . 1c in active galaxies i Pounds & Reeves! 
120091 : iReeves et alj|2009l ; iTombesi et alj|2010ah . Indeed, the out- 
flows required in order to inflate the shocked regions of the ob- 
served nebulae may need to have mildly relativistic o utflow ve- 
locities JPakull & Mirionil2002llAbolmasov et al.l2007h. si milar to 
those seen in the Galactic system SS 433 (see lFabrikal2004l and ref- 
erences therein). Outflows with extreme velocities therefore remain 
an intriguing possibility for ULXs, particularly if IGRJ17091- 
3624 turns out to be radiating at super-Eddington luminosities. 

However, the most obvious possibility, assuming that these 
sources do launch strong, super-Eddington outflows along our line 
of sight, is that in contrast to Galactic BHBs, the ionisation state 
of the material is such that it is not readily observable via iron ab- 
sorption. The ionisation state of a photoionised plasma is typically 
quantified as £ = L/nR 2 , where L is the ionising luminosity from 
the source at a distance R and n is the density of the plasma. The 
outflows in GRS 1915+105 typically display ionisation parameters 
of log£ ~ 4. We know how L changes between GRS 1915+105 
and these bright ULXs, and therefore how the ionisation should 
change for a constant outflow density and distance, but the combi- 
nation of nR 2 is not easy to constrain and will determine how the 
ionisation truly evolves. 

Observationally, a lower ionisation outflow does not seem par- 
ticularly likely in these cases as a massive (high column) outflow 
with relatively low ionisation should have very obvious effects else- 
where in the spectrum, around the oxygen transitions (~0.5-1.0 
keV). Indeed, for brief illustration, we insert a moderately ionised 
absorber using the XSTAR photoionisation code into the model 
used for Holmberg IX X- 1 , fixing the ionisation at log £ = 2 and 
requiring that the redshift of the absorber is either consistent or 
blueshifted with respect to the redshift of the host galaxy. Apply- 
ing this absorber to the full continuum model, we find that the 
column density of any moderately ionised material must be Nn 
< 6 x 10 20 atom cm (note that the model used has a lower limit 
of 10 20 atom cm~ 2 to the allowed column densities), compared to 
the columns of ~ 10 23 atom cm -2 seen in the sub-Eddington out- 
flows of GRS 1915+105 and GRO J1655-40. 

It is therefore probably more likely, based on these observa- 
tional constraints, that the ionisation state of any line-of-sight wind 
is such that the bulk of the iron is close to being fully ionised, 
and hence is prevented from imprinting any strong features on 
the spectrum. Inde ed, continuing their detailed study of the wind 
in GRS 1915+105. lNeilsen et alj 42012b present tentative evidence 
that the wind in this source may start becoming overionised at the 
highest luminosities. However, we note that at the ionisation states 
typically displayed by the wind in GRS 1915+105, increasing £ by 
a factor of ~10, in line with the differences in luminosity between 
this source and the ULXs considered here, only decreases the EW 
of the iron features that would be observed by a factor of ~2 (for 
a constant column density), so even accounting for the sub-solar 
metallicity of the host galaxy, this scenario would require that the 
ionisation increase by a substantially larger factor than could be 
accounted for simply by the difference in luminosity. 

Such a large increase in ionisation could potentially occur 
through either a simultaneous decrease in the distance between the 
ionising source and the outflow material, i.e. the launch radius of 



the wind decreasing, or through a simultaneous decrease in outflow 
density. Given that we expect more material to be expelled at higher 
Eddington ratios, we consider the latter evolution to be unlikely. 
However, the former would be contrary to the larger radius invoked 
in the toy model for hi ghly super-Eddington outflows proposed by 
lMiddletone7ai] ( l2rjT l]), and the more general theoretical prediction 
that the radius from which outflows are launch ed by radiation pres- 
sure i ncreases with increasing M ut/ M acc dShakura & Sunvaevl 
jl973j). Therefore, such a large increase in ionisation appears dif- 
ficult to justify on physical grounds. Furthermore, we note that if 
IGR J17091-3624 does turn out to be a strongly super-Eddington 
source, this overionised scenario would be ruled out on observa- 
tional grounds. 

Individually or in combination, all of these effects can increase 
the difficulty of detecting strong line-of-sight outflows via iron ab- 
sorption, even if such outflows are present. Nonetheless, the limits 
that can currently be placed on any persistent iron absorption fea- 
tures in the spectra of Holmberg IX X- 1 and NGC 1 3 1 3 X- 1 , when 
considered in the context of the behaviour of sub-Eddington BHBs, 
are clearly already beginning to tell us something about their accre- 
tion/outflow processes. Any interpretation invoking highly super- 
Eddington accretion with outflows along our line of sight for these 
sources must be able to explain the surprising lack of strong iron 
features in comparison to sub-Eddington BHBs given the expected 
increase in outflow rates. However, given the apparent difficulty 
in explaining the increase in outflow ionisation required for the 
geometry in which the outflows cross our line of sight to remain 
plausible, the lack of strong absorption features appears to im- 
ply that either super-Eddington outflows are not present in these 
sources, which if true would strongly argue agaist the highly super- 
Eddington interpretation, or that if they are present they must be 
directed away from our line of sight to the central accretion flow. 
If weaker outflows are detected in the future with more sensitive 
observations of these sources, this would suggest that the latter sce- 
nario is not the case. 



4.2 Emission Lines 

Until this point, we have been discussing the potential presence of 
the strong outflows associated with highly super-Eddington accre- 
tion specifically assuming the ejected material impinges our line of 
sight to the central accretion flow. We now consider the possibility 
that strong super-Eddington outflows are present in Holmberg IX 
X-l and NGC 1313 X-l but, despite the large solid angles such 
outflows are expected to cover, are directed away from our line of 
sight. In this scenario, one might expect to see evidence of emission 
features as the outflowing material absorbs and reprocesses X-rays 
from the central source. Narrow iron emission features, particularly 
from neutral iron, are ubiquitously observed in the X-ray spectra 
of Galactic HMXBs, as the strong stellar winds from the massive 
companions reproce ss the X-rays from the central black hole (see 
iTorreion et al.ll2010r) . As with HMXBs, the stellar companions of 
ULXs are also expected to be fairly young, relati vely massive stars 
(confirmed in the case of Holmberg IX X-l bv lGrise et af]|201ll) 
which should naturally eject strong stellar winds regardless of the 
nature of the compact object powering the ULX. 

Again considering the immediate Fe K band, we find very sim- 
ilar constraints to those placed on absorption features; any narrow 
emission features present must have equivalent widths less than 
~30and ~50 eV for Holmberg IX X- 1 and NGC 1313 X-l respec- 
tively. Furthermore, emission features should not display the sig- 
nificant velocity shifts possible for absorption features, so in prin- 
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ciple we can be more precise about the exact energies considered 
(note that any additional peculiar velocities of these sources within 
their host galaxies could lead to the rest frame energy scale adopted 
here being systematically offset, but such effects are unlikely to be 
large). Assuming the redshifts quoted in Table[T]are appropriate for 
these sources as well as their host galaxies, the 99 per cent limits on 
emission features obtained for the transitions highlighted in Figs. [2] 
and [3] (neutral, helium-like and hydrogen-like iron) are EW < 23, 
25 and 12eV for HolmberglX X-l, and EW < 41, 37 and 11 eV 
forNGC 1313 X-l respectively. 

From these equivalent width limits, we can obtain upper lim- 
its for t he reprocessin g colu mn density using the relation de- 
rived by iKallman et ail J2004b: EW = ANh, w here A is some 
constant of proportionality. IKallman et al. I J2004h showed that for 
neutral iron, A neutra i ~ 3 (for a solar iron abundance, equiv- 
alent widths in eV and column densities in 10 22 atom cm -2 ), 
w hich has been observa tionally confirmed for Galactic HMXBs 
bv lTorreion et al. Using this relation, (although correcting 

for the ~0.5 Zq metallicities of the host galaxies) we find that 

neutral 

< 14 and 24 x 10 22 atom cm for HolmberglX X- 
1 and NGC 1313 X-l respectively, placing both at the lower end 
of the distribution observed in Galactic HMXBs, which ranges 
from~10 22 - 10 24 atom cm , with the majority disp laying Nn 
> 10 23 atom cm -2 ; see Fig. 5 in lTorreion et aljfeoid . However, 
that work also reported an anti-correlation between the neutral iron 
equivalent width and HMXB luminosity, which they interpret as be- 
ing due to an increase in ionisation of the stellar winds in which the 
accreting black holes are embedded. The limits obtained here for 
neutral iron are consistent with this picture. Any super-Eddington 
outfows present can not provide any substantial neutral column in 
addition to the stellar wind expected from the companion, but given 
the high ionisation of sub-Eddington outflows in BHBs, we would 
naturally expect any super-Eddington outflows to be at least mod- 
erately ionised. 



If the outflows have a similar ionisation to the outflows ob- 
served in GRS 1915+105 (log£ ~ 4), iron will be primarily 
in its hydrogen-like ionisation state. Re-calculating the constant 
of proportionality A for Fe XXVI, taking the threshold absorp- 
tion cross-section of hyd rogenic iron to be ~9.5 xl0~ 21 cm 2 
dVerner & Yakovlevll 19951) , the iron abundance to be ~ 1.6 xlO 



(e.g. I Wilms et al ] |200Ct Z ~ 0.5 Zq), and the flourescent yield to 
be ~0.5, we find that iV H ,iogf~4 < 9 and 8 x 10 22 atom cm -2 for 
Holmberg IX X- 1 and NGC 1313 X-l respectively. Therefore, there 
does not appear to be any substantial column of highly ionised re- 
processing material either. 

However, we return to the issue that the expected evolution of 
the ionisation state of the wind with increasing Eddington ratio is 
not well known, owing to the uncertain evolution of the combina- 
tion nR 2 . As discussed previously, an increase in ionisation will 
increase the column density inferred from the limit on the emission 
from Fe XXVI as a larger fraction of the iron nuclei become fully 
ionised. However, in comparison to the line-of-sight (absorption) 
case, it is more observationally plausible for the ionisation to have 
decreased in this scenario. Indeed, as the Eddington ratio increases, 
standard theory predicts that the launch radius of the wind should 
also increase, and it is not unreasonable to imagine the outflow den- 
sity might also increase; both of these effects could serve to coun- 
teract the increase in luminosity and result in an overall decrease in 
ionisation. If super-Eddington outflows are sufficiently less ionised 
than sub-Eddington outflows, they could enter a regime in which 
resonant trapping and Auger line destruction is efficient, which 



would result in only very weak emission lines (see e.g. iRoss et all 
[l996h . 

A further caveat to our consideration of emission lines is that 
the relat ion between line equi valent width and column density de- 
rived by |KallmanetalJ ( l2004h assumes a spherical distribution for 
the reprocessing medium. For a constant column, a non-spherical 
outflow distribution (i.e. covering a solid angle < 47r) would 
reduce the EW expected for the reprocessed emission lines. Al- 
though a spherical distribution is certainly appropriate for the con- 
sideration of reprocessing by any stellar winds from the compan- 
ion, we now consider whether it is also a reasonable assumption for 
strong super-Eddington outflows. It is widely expected that as the 
Eddington ratio i ncreases, the scale height of the wi nd should also 
increase (see e.g. King 20091 iDotan & ShavhfeOl lh . so as the Ed- 
dington ratio increases, the assumption of a spheri cal wind distribu- 
tion be comes more and more appropriate. Indeed, fbotan & Shavivl 
d20 1 lh find that the winds launched by super-critical accretion be- 
come spherical beyond L/Le — 5, although this is somewhat de- 
pendent on the assumed model parameters, e.g. the disc viscosity 
and the opacity parameterisation. Nevertheless, although by virtue 
of the lack of strong absorption features any strong outflows cannot 
be truly spherical in these cases, the high Eddington ratios (~10) 
inferred for black hole masses of 10 Mq might well result in some- 
thing approaching a spherical outflow geometry, so the estimated 
limits on the columns obtained might not be unreasonable. We 
stress that while reducing the solid angle covered by the outflow 
is a legitimate way in which to potentially increase the limits on 
the column density of the outflow, doing so would also imply a 
reduction in the true Eddington ratio. 



4.2.1 Broader Lines 

In this work, we have focused primarily on narrow iron features, 
given the observed absorption features from sub-Eddington out- 
flows and the emission features from Galactic HMXBs. How- 
ever, the emission lines expected from an optically thick, diverg- 
ing outflow may be broadened by both Doppler shifts and elec- 
tron _sj^attermg_ 2 _pju'ticjiW velocity is large (see 
e.g. iTitarchuk et alj|2009l[slm et alj|20ld) . Therefore, we consider 
briefly the limits that can currently be placed on emission lines as a 
function of the line width. For illustrative purposes, we demonstrate 
in Fig.|4]how the EW limits evolve over a range of assumed line 
widths for the K a transitions of Fe I, XXV and XXVI for Holm- 
berg IX X-l. The limits naturally become weaker as the assumed 
width increases, therefore stronger iron emission features could yet 
be present if they are substantially broadened by the wind in which 
they arise. 

Unfortunately, the combination of the uncertainty in the ion- 
isation state, solid angle and the outflow velocity of the wind, all 
of which contribute towards determining the width of any emission 
line produced by the outflow, means it is difficult to make any clear 
predictions about the expected line widths, and hence provide a de- 
tailed consideration of the limits on any such emission lines. We 
therefore conclude that at this stage it is still too difficult to place 
strong constraints on the presence of super-Eddington outflows that 
do not cover our line-of-sight just through consideration of the ex- 
pected iron emission. 
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Figure 4. As in Fig. [2] but here we plot the evolution in the EW limits with the assumed line width. For illustrative purposes, we only show evolution of 
the limits obtained at the energies of the neutral, helium-like and hydrogen-like iron transitions for Holmberg IX X- 1 . Larger equivalent widths are naturally 
allowed as the line width increases, as by broadening the line the overall line flux can be further increased without substantially perturbing the fit to the data. 



5 CONCLUSIONS 

Making use of all the good quality XMM-Newton data publically 
available for two bright ULXs, Holmberg IX X-l and NGC 1313 
X-l, we have searched for (persistent) discrete atomic features in 
their high energy spectra that could be associated with either iron 
emission or absorption, and provide observational evidence for the 
massive outflows predicted if these sources are accreting at substan- 
tially super-Eddington rates. We do not find any statistically com- 
pelling evidence for any features, either in absorption or in emis- 
sion. However, we do find that the data currently available for these 
sources is of sufficient quality that we can place fairly stringent 
limits on the strengths of any features intrinsically present these 
spectra that remain undetected. Any features present (absorption or 
emission) in the immediate Fe K energy band (6-7 keV) must have 
equivalent widths weaker than ~30eV for Holmberg IX X-l, and 
weaker than ~50eV for NGC 1313 X-l. 

In comparison to the strongest sub-Eddington outflows ob- 
served in GRS 1915+105, which imprint iron absorption features 
with equivalent widths of ~30 eV, the limits obtained for the ULXs 
considered here appear quite restrictive, particularly when these 
sources must be radiating at ~5-10 times their Eddington limits if 
they host black holes of similar masses, and should therefore be ex- 
pelling at least 5-10 times as much material as GRS 1915+105. The 
difficulty in trying to reconcile these observational constraints with 
the presence of strong line-of-sight outflows leads us to conclude 
that either these sources do not launch such outflows, which would 
strongly argue against a highly super-Eddington interpretation, or 
that they much be launched away from our viewing angle. Addi- 
tional deep observations of bright ULXs with current instrumenta- 
tion, and in particular with the micro-calorimeter due for launch on 
board Astro-H in the near future, will be essential for detecting any 
weak iron features in such sources, and hence in placing stronger 
constraints on the nature of any outflows present, and the ultimate 
nature of these sources. 
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APPENDIX A: INDIVIDUAL OBSERVATION DETAILS 

Here, we examine briefly the individual observations of Holm- 
berg IX X-l and NGC 1313 X-l considered in this work to check 
whether it is reasonable to combine the data obtained into a single 
spectrum in each case. For each source, we model the spectra ob- 
tained from the individual observations simultaneously, adopting 
the same continuum model outlined in the main body of the pa- 
per (see section [3~Tt . The key continuum parameters are all free to 
vary between the observations, although we require that the neutral 
absorption is the same for each. Owing to the highly varied data 
quality it was not always possible to reliably constrain both the op- 
tical depth and electron temperature of the Comptonising region 
in NGC 1313 X-l, as these parameters often display some degen- 
eracy. In these cases, we fixed the electron temperature at 2 keV, 
similar to the values obtained for the other observations. The results 
obtained are presented in Table lAll Although there are some minor 
differences in the parameter values obtained, they are all reasonably 
similar for both sources. 

In addition, to provide a visual comparison of the spectra from 
the various observations, we plot the count spectra of the different 
observations for each source in Fig. I A 1 1 This confirms the results 
obtained with the quantitative analysis presented in Table lAll there 
are some slight differences between the spectra of the individual 
observations in terms of their high energy spectral slope/curvature, 
more so for NGC 1313 X-l than HolmberglX X-l, for which the 
high energy spectra are practically identical, but ultimately the 
spectra are all fairly similar for both sources. Hence we conclude 
that it is reasonable to combine them into a single, averaged spec- 
trum for each source in order to improve the high energy statistics 
and further facilitate our search for narrow iron features. 

As a final test, for each of the models obtained for the ob- 
servations of NGC 1313 X-l (see Table lAll we inserted a narrow 
(a = 10 eV) Gaussian absorption line at 6.67 keV with an equiv- 
alent width of -50 eV. Data were then simulated from each model 
using the EPIC-pn response, keeping the respective exposures of 
the simulated spectra the same as the good exposure times for the 
observations on which they are based (see Table[TJ. We then com- 
bined the simulated datasets in the same manner as the real data, 
and examined the characteristics of the absorption feature in the 
combined simulated dataset. The narrow feature is well detected, 
providing an improvement of A% 2 = 20 for 2 additional degrees 
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Table Al. Key parameters obtained for all the individual observations of HolmberglX X-l and NGC 1313 X-l with the phenomenological ultraluminous state 
continuum model (see text). 



Source 


Obs. 


JV H 


kT in 


kT c 


T 


Xl (X 2 /d-o.f.) 






(10 21 cm" 2 ) 


(keV) 


(keV) 







HolmberglX X-l 



NGC 1313 X-l 



0112521001 
0200980101 
0200980101 

0106860101 
0150280601 
0150181101 
0205230201 
0205230301 
0205230601 
0405090101 



1.5 ± 0.1 



2.6 ±0.1 



0.24 ±0.04 


2.6+"-5 

— U.4 


7.5 ± 1.2 


1.02(3162/3102) 


0.22 ± 0.05 


2 9+ 1 ' 9 






0.24 ±0.01 


2.3 ±0.1 


9.2 ± 0.4 




0.21 ± 0.02 


9 n+0.4 
z - -0.2 


8.2 ± 0.9 


1.05 (3968/3769) 


0.30 ± 0.06 
S5+ 019 


2* 
1 6+ 1 ' 2 


7 R+ 1 ' 1 

> 5.3 




0.25 ± 0.05 


2* 






a28 ^o.o<3 

n 9 o+0.03 
u ' zo -0.02 


1 6 +0 ' 2 
i - D -0.3 


8 4+ 3 ' 4 
q o+ - 6 




0.22 ± 0.01 


n 9 + 0.1 
z - -0.2 


8.5 ± 0.5 






I 10 1 10 



Energy (keV) Energy (keV) 

Figure Al. Count spectra for the individual observations of Holmberg IX X-l (left panel) and NGC 1313 X-l (right panel), and the data/model ratios for the 
best fit continuum models presented in Table lAll For clarity, only the EPIC-pn data are shown, and the data are rebinned for display purposes only. 



of freedom, and we find a line energy of E = 6.65 ± 0.04 keV, 
and an equivalent width of EW = —52 ± 19 eV (90 per cent con- 
fidence). Considering instead the simulation of the long observa- 
tion by itself, the absorption feature only offers an improvement of 
A\ 2 = 14 for 2 extra degrees of freedom. Therefore, even given 
the mild variance in the high energy spectrum of NGC 1313 X-l, 
combining the individual spectra does aid the detection of persis- 
tent narrow features. 



© 2010 RAS, MNRAS OOO.rflfiTl 



